rj and r]' mesons in the Dyson-Schwinger approach using a generalization of the 

Witten-Veneziano relation 
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The description of the ry and ry' mesons in the Dyson-Schwinger approach has relied on the Witten- 
Veneziano relation. The present paper explores the consequences of using instead its generalization 
recently proposed by Shore. On the examples of three different model interactions, we find that 
irrespective of the concrete model dynamics, our Dyson-Schwinger approach is phenomenologically 
more successful in conjunction with the standard Witten-Veneziano relation than with the proposed 
generalization valid in all orders in the 1/Nc expansion. 

PACS numbers: ll.lO.St, 12.38.-t, 14.40.Aq, 12.38.Lg 



I. INTRODUCTION 

The Dyson-Schwinger (DS) approach [H, i, i, i, i, Q 
to QCD and its modeling is the chirally well-behaved 
bound-state approach. Thus, it is the most suitable 
bound-state approach to treat the hght pseudoscalar 
mesons (those composed of the u,d and s quarks), for 
which dynamical chiral symmetry breaking (DChSB) is 
essential. One solves the ("gap") DS equations (DSEs) 
= {S^^"^)^^ ~ for tire dynamically dressed 
quark propagators Sg, where q is the quark flavor {q = 
u,d,s, ...), Sg '^'^ is the free quark propagator, and is 
the quark self-energy. These dressed quark propagator 
solutions are then employed in Bethe-Salpeter equations 
(BSEs) for the boimd-state vertex Tqq> of the meson com- 
posed of the quark of the flavor q and antiquark of the 
flavor q': 



[^qq']ef — j [SqT qqi Sqi]gh 



(1) 



where e,f,g,h schematically represent spinor, color and 
flavor indices, integration is meant over loop momenta, 
and K is the interaction kernel. Solving Eq. ^ for Tqq' 
also yields Mqq> , the mass eigenvalue of the qq' meson. 

To obtain the chiral behavior as in QCD, DS and 
BS equations must be solved in a consistent approxima- 
tion. The rainbow-ladder approximation (RLA), where 
DChSB is well-understood, is still the most usual ap- 
proximation in phenomenological applications. This also 
entails that in both DSE and BSE ([T]) we employ the 
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same effective interaction kernel. 



A'' 



[Kiktf = ic/D'/iikUiYrU [yTl./ , (2) 

so that the quark self energy in the gap DSE is 

^'^(p^ = - l-0j49'D';i{kU^rSqie)^r ■ (3) 

In Eqs. ^ and ([3]), £'|J^(fc)off is an effective gluon propa- 
gator. For example, for renormalization- grou p improved 
(RGI) interactions (e.g., in Refs. BBSMflii)- it has 
the form 



(4) 



where D°(^(p)free is the free gluon propagator, and 
aeff(fc^) is an effective running coupling. For large space- 
like momenta {k^ ^ 1 GeV^), aoff(fc^) approaches the 
perturbative QCD running coupling as{k'^) known from 
the QCD renormalization group analysis, although it 
must be modeled at low momenta. 

Concretely, in the present paper we recall and utilize 
the results obtained i) in Refs. 0, [1] by using the RGI 
of Jain and Munczck 0, ii) in Ref. by using the 
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RGI gluon condensate- induced interaction and Hi) 
in Refs. [H, [I!] by using the separable interaction [TJl- 
In any case, such effective interactions must be modeled 
at least in the low-energy, nonperturbative regime in or- 
der to be phenomenologically successful - which above 
all means to be sufficiently strong in the low-momentum 
domain to yield DChSB. In the chiral limit (and close 
to it), light pseudoscalar (P) meson qq bound states 
{P = TT*''^, 77) then simultaneously manifest them- 

selves also as (gitasi-) Goldstone bosons of DChSB. This 
enables one to work with the mesons as explicit qq bound 
states, while reproducing the results of the Abelian axial 
anomaly for the light pseudoscalars, i.e., the amplitudes 
for P 77 and 7* P'^P^P~. This is unique among 



2 



the bound state approaches - e.g., see Refs. [sl. [isl. [TgI. [TtI 
and references therein. Nevertheless, one keeps the ad- 
vantage of bound-state approaches that from the qq sub- 
structure one can calculate many important quantities 
(such as the pion, kaon and ss pseudoscalar decay con- 
stants: /tt, Jk and fss) which are just parameters in most 
of other chiral approaches to the light-quark sector. The 
treatment 7, 8, 10,^8,] of the r/-?/ complex is remarkable 
in that it is very successful in spite of the limitations 
of RLA. (Very recently, during the work on the present 
paper, the first and still simplified DS treatments of rj 
and 77' beyond RLA appeared [l^, H^]. However, RLA 
treatments will probably long retain their usefulness in 
applications where simple modeling is desirable, as in 
the calculationally demanding finitc-tcmperature calcu- 
lations [13.) The RLA treatments of the 77-77' complex 
at first determined 0, H, [ill ^^^^ anomalous 7^0 mass pa- 
rameter by fitting the empirical ry and 77' masses. More 
recently, the treatment was improved by avoiding this 
fitting while retaining the phenomenologically successful 
description • Namely, the anomalous 770 mass was 

no longer a free parameter but determined from the lat- 
tice results (on QCD topological susceptibility) through 
the Witten-Veneziano (WV) relation [2l|, [23|. However, 
Shore achieved [l^, [l^l what can be considered as a gen- 
eralization of the WV relation, and the purpose of the 
present paper is exploring the usage of this generaliza- 
tion in the DS context. 

The paper is organized as follows: in the next sec- 
tion, we recapitulate the procedures and results of our 
previous treatments [1, [To, relying on the WV rela- 
tion ([T7)) . and present in Table I also their extension to 
the scheme of the four decay constants (and two mixing 
angles) of -q and 77'. In Section Hill w e expose the usage 
of the pertinent Shore's equations |23l. [2^ in the context 
of DS approach. The last section concludes after giving 
the results of solving the pertinent equations. 



II. T]-T)' MASS MATRIX FROM 
WITTEN-VENEZIANO RELATION 

All qq' model masses Mqqi {q, q' = it, d, s) used in the 
present paper, and corresponding qq' bound-state ampli- 
tudes, were obtained in Refs. |3, H, [l^j El, dj HI] in 
RLA, i.e., with an interaction kernel which (irrespective 
of how one models the dynamics) cannot possibly cap- 
ture the effects of the non-Abelian, gluon axial anomaly. 
Thus, when we form the 77-77' mass matrix 



Mis Af|o 
Mis 



(5) 



in this case in the octet-singlet basis 778-770 of the (broken) 
flavor-SU(3) states of isospin zero, 

r\% — ^={uu dd ~ 2ss) ^ ?]o ^ —={uu + dd + ss), (6) 
v6 v3 



this matrix ([5]), consisting of our calculated qq masses. 
Mi, ^ ivslM^M = I {M% + \mI^) , (7) 



M|o = (778|MSaI%> = Mis = ^(Af|^_^_A4) < Q, (8) 



M. 



2 _ 
00 — 



(9) 



is purely non-anomalous (NA), vanishing in the chi- 
ral limit. In the isospin limit, to which we adhere 
throughout, the pion is strictly decoupled from the gluon 
anomaly and M-^u = M^j is exactly our model pion mass 
M^. Also the unphysical ss quasi-Goldstone's mass Msg 
results from RLA BSE and does not include the con- 
tribution from the gluon anomaly. This is consistent 
with the fact that due to the Dashcn-GcU-Mann-Oakes- 
Renner (DGMOR) relation, it is in a good approximation 
BiEIil given by 



Mt 



2M 



K 



Mi 



(10) 



i.e., by the kaon and pion masses protected from the 
anomaly by strangeness and/or isospin. 

In our previous DS studies 0, i, 0, [H, [H, HI, to 
which we refer for all model details, the phenomenology 
of the non-anomalous sector was successfully reproduced, 
e-g-, f-n, Jk, as well as the empirical masses M^r and Mk 
(see the upper part of Table |T|, yielding a strongly non- 
diagonal M^^ ([5]) . Its diagonalization leads to the cigen- 
states known as the nonstrange-strangc (NS-S) basis. 



77ns = -^iuu + dd) 



ss . 



(11) 



and to M^j^ = diag[Af|, M|j]. In contrast to these mass- 
squared eigenvalues, the experimental masses are such 
that [Ml)exp < (M2)e^p, and 77' is too heavy, {M,y)exp = 
958 MeV, to be considered even as the ss quasi-Goldstone 
boson. This is the well-known Ua{^) problem, resolved 
by the fact that the complete rj-rj' mass matrix M^ must 
contain the anomalous {A) part M\. That is, M^ = 



M2 



Ml. 



However, M^ is inaccessible to RLA which yields our 
Goldstone pseudoscalars. In Refs. 0, B [H, 111 , M| 
was extracted from lattice data through the WV relation 
[the second equality in Eq. p7|) ]. The purpose of the 
present paper, instead, is to approach 77 and 77' through 
Shore's [23, [2J recent generalization of that relation. 

Before that, however, we review the usage of the 
WV relation in Refs. 0, i, [13, [11 [11 ■ The expan- 
sion in the large number of colors, Nc, indicates that the 
leading approximation in that expansion describes the 
bulk of main features of QCD. The gluon anomaly is sup- 
pressed as 1/A^c and can be viewed as a perturbation in 
the large expansion. In the SU(3) limit [compare Eqs. 
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((T2|) and ([T3)) ] , it is coupled otiZj/ to the singlet combina- 
tion 770 (O; only the tjq mass receives, from the gluon 
anomaly, a contribution which, unlike quasi-Goldstonc 
masses Mggi's comprising Af^^, does not vanish in the 
chiral limit. As discussed in Refs. 0, the present 

bound-state context it is thus meaningful to include the 
effect of the gluon anomaly just on the level of a mass 
shift for the 770 as the lowest-order effect, and retain the 
qq bound-state amplitudes and the corresponding mass 
eigenvalues M^^ as calculated by solving DSEs and BSEs 
with kernels in RLA. 

References 0, i, [13, d, [H thus break the C/a(1) 
symmetry, and avoid the J7a(1) problem, by shifting the 
r)o (squared) mass by an amount denoted by 3/3 (in the 
notation of Refs. [1, [l^l)- The complete mass matrix 
Af^ ~ M^p^ + M\ then contains the anomalous part 



Ml = diag[0,3/3] 



(12) 



where the anomalous 770 mass shift 3/3 is related to 
the topological susceptibility of the vacuum, but in the 
present approach must be treated as a parameter to be 
determined outside of our RLA model, i.e., fixed by phe- 
nomenology or taken from the lattice calculations [26| . 
(The possibility of employing an additional microscopic 
model for the gluon anomaly contribution, such as the 
one of Ref. is presently not considered.) 

The SU(3) flavor symmetry breaking and its inter- 
play with the gluon anomaly modifies M\ (fT2|) to 



Ml 



f (1 



3 



X){2 + X) 



{l-X){2 + X) 
\{2 + Xf 



(13) 

where X is the fiavor symmetry breaking parameter. It is 
most often estimated as X = Jt^I fss ^ 0.7 — 0.8 (see^e.g., 
Refs. H, [13, Hill, although there are some other of 
course related, estimates of X). Presently we also adopt 
X = f-n / fssi which means that A" is a calculated quantity 
in our approach. The employed models achieved good 
agreement with phenomenology 0, 0, [13, ES] , e.g., fitted 
the experimental value of Al"^ + M^, for (3 around 0.26 - 

0.28 GeV^. The anomaly contribution M\ then brings 
the complete rather close to a diagonal form for all 
considered models 0, H, [13, [II]; that is, to diagonalize 
M^, only a relatively smaU rotation {\e\ - 13° ± 2°) of 
the Tj^-r^Q basis states. 



T] = cos U rjs ~ sm rjo 



77' = sin 6 r/8 4- cos 9 r/o , (14) 



is needed to align them with the mass eigcnstates, i.e., 
with the physical 77 and 77'. In contrast to this, the 7;- 7;' 
mass matrix in the NS-S basis (|lip . 



ms nsvNs 

msvs vs 
Ml + 2(3 V2PX 

V2f3X M% + f3X^ 



(15) 



Af 2 
A4 



is then strongly off-diagonal. The indicated diagonaliza- 
tion, given by 

77 = cos 77ns - sin (/> 77s , 77' = sin (f> t^ns + cos (f> ijs , 

(16) 

is thus achieved for a large NS-S state-mixing angle </) ~ 
42° ±2°. Of course, this is again in agreement with 
phenomcnological requirements [l3 | , since <f> is fixed 
to the ris-Tjo state-mixing angle 6 by the relation (f> ~ 
6 + arctan \/2 = 6 + 54.74°. The masses are 



M^ = - 
V 9 



AC = - 
2 



A/i 



A/i 



A/i 



(ML 



- Af 2 )2 + 8/32^^ 



The invariant trace of the mass matrix (|15p . together 
with Eq. (dni), gives the first equality in 



13 {2 + X'^) = A/2 + A/2 



2Ml 



K 



6 

72 

J TT 



(17) 



The second equality is the Witten-Veneziano (WV) re- 
lation [U, [l^l between the 77, 77' and kaon masses and 
Xym, the topological susceptibility of the pure gauge, 
Yang-Mills theory. Thus, [3 does not need to be a free 
parameter, but can be determined from lattice results on 
Xym, so that no fitting parameters are introduced. For 
the three models [3, [U, [3 utilized in our treatments 
0, [1, [13, [13] of 77 and 77', the bare quark mass parameters 
and the interaction parameters were fixed already in the 
non-anomalous sector, by requiring the good pion and 
kaon phenomenology. (See the tt and K masses and de- 
cay constants in the uppermost part of Table HI) Then, 
following Refs. [l3, [ll] in adopting the central value 
of the weighted average of the recent lattice results on 
Yang-Mills topological susceptibility [33, [Ml, [13] , 



XYM = (175.7 ± 1.5 MeV)^ 



(18) 



we have obtained the good descriptions of the 77-77' phe- 
nomenology 0, [3, [13, [13] , exemplified by the first three 
columns (one for each DS models used) of the middle 
part of Table U giving the predictions for the rj and rj' 
masses and for the NS-S mixing angle (/). 

The lowest part of the table, below the second hori- 
zontal dividing line, contains the results on the quantities 
(6*0, ^8, etc.) defined in the scheme with four 77 and 77' de- 
cay constants and two mixing angles, introduced and ex- 
plained in the following Section lTlll TableUalso compares 
these results of ours (in the first three columns) with 
the corresponding results of Shore's approach [13, [13, 
which the experimental values of the meson masses A/^, 
Mk, Mrj, and A/^/, as well as the decay constants and 
Jk (in contrast to our qq bound-state model predictions 
for these quantities) are used as inputs enabling the cal- 
culation of various decay constants in the 77-77' complex 
and the two mixing angles a-nd 63 (corresponding to 
(j> — 38.24° in our approach). 
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III. USAGE OF SHORE'S EQUATIONS IN DS 
APPROACH 

The WV relation was derived in the lowest-order ap- 
proximation in the lar ge N r. expansion. However, consid- 
erations by Shore [23l |24| contain what amounts to the 
generalization of the WV relation, which is valid to all 
orders in 1 /Nc ■ Among the relations he derived through 
the inclusion of the gluon anomaly in DGMOR relations, 
the following are pertinent for the present paper: 

(/°,)^A/,^^ + iO'M'^ = lif^M'^ + 2/iMi) + 6A , (19) 
/° /8 + /0/«m2 ^ ^ (/^M^ _ /^M^) , (20) 
{f?r)'M?r + ifnfMf, = - \ [llMl - 4/i Afi) , (21) 



from 


[8J 


[10] 


[13] 


Shore 




Ref. 


& WV 


& WV 


& WV 


[23, 2^ 


Experiment 




137.3 


135.0 


140.0 




(-1 OQ ri\isospin 
V -'-'-'^•^/ aver age 


Mk 


495.7 


494.9 


495.0 




('495 vysospin 

V^'-"-'' ' ^average 




700.7 


722.1 


684.8 






U 


93.1 


92.9 


92.0 




92.4 ±0.3 


fK 


113.4 


111.5 


110.1 




113.0 ± 1.0 


fss 


135.0 


132.9 


119.1 






Mr, 


568.2 


577.1 


542.3 




547.75 ± 0.12 


M„, 


920.4 


932.0 


932.6 




957.78 ±0.14 


4> 


41.42° 


39.56° 


40.75° 


(38.24°) 




Bo 


-2.86° 


-5.12° 


-6.80° 


-12.3° 




es 


-22.59° 


-24.14° 


-20.58° 


-20.1° 




fo 


108.8 


107.9 


101.8 


106.6 




fs 


122.6 


121.1 


110.7 


104.8 




rO 

Jv 


5.4 


9.6 


12.1 


22.8 




Jr,' 


108.7 


107.5 


101.1 


104.2 




Jv 


113.2 


110.5 


103.7 


98.4 




Jv' 


-47.1 


-49.5 


-38.9 


-37.6 





TABLE I: The results of employing the WV relation (ITTll 
in our DS a ppr oach for the three dynamical models used in 
Ref s. jsl. [lol [l3 | . compared with the results of Shore's analy- 
sis [23l . [24 1 and with the experimental results. The first col- 
umn was obtained by the WV-recalculation of the results of 
Ref. which in turn used the Jain-Munczek Ansatz for 

the gluon propagator 0]. Column 2: the results based on 
Ref. [l3|, which used the OPE-inspired, gluon-condensate- 
enhanced gluon propagator [l3| . Column 3: the results based 
on Ref. uMj which utilized the separable Ansatz for the 
dressed gluon propagator [ij. Column 4: The results of 
Shore [23. 12411 . who used the lattice result XYM = (191 MeV)** 
of Ref. [3^, and not the weighted average |(T5}, in contrast 
to us. Column 5: the experimental values. All masses and 
decay constants are in MeV, and angles are in degrees. For 
more details, see text. 



where A is the full QCD topological charge parameter, 
and /°,, /°, /^,, are the four decay constants [H, 113, 
Issj associated with the two isoscalar pseudoscalars 77 and 

The nonperturbative parameter A is related to the 
QCD topological susceptibility, quark condensates and 
quark masses [2^ At large Nc, it should be well- 

approximated by the topological susceptibility, A ^ x- 
More precisely, it reduces to the YM topological suscep- 
tibility in the large Nc limit: A = xym + 0{1/Nc), but at 
present it is not known better than that, as there are still 
no lattice data on this nonperturbative QCD parameter. 
Therefore, in his own phenomenological analysis. Shore 
himself had to approximate A by a value of xym [H, HI] . 
In that sense, because of this crucial assumption based 
on the lowest-order 1/Nc approximation, even his analy- 
sis was not (and, because of the lack of the corresponding 
lattice data, could not be) carried out numerically con- 
sistently in the orders of Nc, even though his formulas 
are valid in all orders in the 1 / Nc expansion. 

While the present bound-state DS approach clearly 
cannot improve on the consistency aspect, it offers the 
possibility of a phenomenological analysis entirely differ- 
ent from Shore's. Namely, in addition to A w Xym, Shore 
used the experimentally known quantities (pion, kaon, ry 
and 77' masses, as well as the pion and kaon decay con- 
stants) as inputs in Eqs. (|19p - ([2T|) to obtain the 77 and 
77' decay constants f^> , f^, f^> , f^- On the other hand, 
the predicting power of our bound-state DS approach is 
much larger: not only are pion and kaon masses and de- 
cay constants calculated quantities, predicted from the 
qq substructure, but once we formulate the incorpora- 
tion of Shore's generalization within the bound-state DS 
approach, it will become obvious that also these four 77 
and 77' decay constants and their masses and Mjfi 
come out as pure predictions. Such a phenomenological 
analysis, complementary to Shores, motivates us to for- 
mulate and perform the treatment based on Shore's gen- 
eralization, instead of the original WV relation (or fitting 
the anomalous 770 mass shift) as in our earlier references 

Adding Eqs. (fT9| and (f2T|) . one gets the relation 

(/°)^A4 + (/°)^M,^ + (4)^M,^ 

± (/«)'m2 -2/iMi = 6A (22) 

which is the analogue of the standard WV formula (I17p . 
to which it reduces in the large Nc limit where A Xym , 
the , /,p fx — > f-w limit, and the limit of vanishing sub- 
dominant decay constants (since rj and 77' are dominantly 
778 and 770, respectively), i.e., f^,f^, 0. However, we 
will need to use not just this single equation, but the 
three equations (fT9| - ([2T|) from Shore's generalization. 

These four 77 and rf decay constants are often pa- 
rameterized in terms of two decay constants, /g and /o, 
and two mixing angles, 0s and 

= cos 08 /8 , /° = - Sin 00 fo , (23) 
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f^, = sin 6*8 fa 



J if 



cos 6*0 /o 



(24) 



This is the so-caUed two-angle mixing scheme, which 
shows exphcitly that it is inconsistent to assume that the 
mixing of the decay constants foUows the pat tern (I14p of 
the mixing of the states rj^ and rjQ [H, [29l. l33l. [3^ . l35l. l36l. 

The advantage of our model is that, as wc shall see, 
we are able to calculate the /§ and fo parts of the phys- 
ical decay constants (|23|l -p4 |) from the qq substructure. 
However, we cannot keep the full generality of Shore's 
approach, which allows for the mixing with the gluonic 
pseudoscalar operators, and therefore employs the defi- 
nition [2^, [i^l of the decay constants which, in general, 
due to the gluonic contribution, differs from the follow- 
ing standard definition through the matrix elements of 
the axial currents A°-^{x): 



{0\A'^'^ix)\P{p))=tf^p^- 



-ip-x 



= 8,0; P = v,7f . 

_ _ (25) 
Nevertheless, Shore's definition [23, [2^ coincides with 
the above standard one in the non-singlet channel, where 
there cannot be any admixture of the pseudoscalar glu- 
onic component. Similarly, since our BS solutions (from 
Refs. 0,1 Em) are the pure qq states, without any 
gluonic components, using Shore's definition would not 
help us calculate the gluon anomaly influence on the de- 
cay constants. We thus employ the standard definitions 
also used by, e.g., Gasser, Leutwyler, and Kaiser 
[33 . [34I. [33 . a s well as by Fcldmann, Kroll, and Stech 

(FKS) [Mli,!!^. 

In Eqs. ((23|)-([24|), the angles are chosen so ^ that 
9s = 00 ^ 9 = in the limit of the exact SU(3) fla- 
vor symmetry, since only then there are just two de- 
cay constants, purely octet = /g and purely singlet 
fry = /o, while the off-diagonal decay constants vanish, 
/° = = /^/, in this limit. Otherwise, all four decay con- 
stants ([25]) are different from zero due to the breaking of 
the SU(3) flavor symmetry, since this leads to 6* 7^ and 
gives both rj and rj' the both components 778 and 770. [In 
the parameterization ([23|) - p4| . the angles 6*8 and 9q dif- 
fer from 9 since also {0\Al\r]o) ^0 ^ {0\A°^\tjs) ■] Thus, 
although not r]s but tjq couples to the gluon anomaly, the 
octet-chanel constants ff^ and f^, are influenced by the 
gluon anomaly through its interplay with the SU(3) fla- 
vor symmetry breaking [similarly to the anomalous mass 
matrix ([T3|) having nonvanishing 88, 08 and 80 elements 
when X ^1]. 

Equivalently to f°,,f^,f°, and f^y, defined by 
Eq. (|25p . one has four related but different constants 
/^^, /^^, and /^/, if instead of octet and singlet axial 
currents (a = 8, 0) in Eq. p5|) one uses the nonstrange- 
strange axial currents (a — NS, S) 



W = -^A'^{x) + ^A'^ix) = ±s{x)ri5s{x) . 

(27) 



The relation between the two equivalent sets is thus 



J rj J rj 
fNS fS 
J rj' J Tj' 



J Tj J rj 
J rj' J rj 



1 

V3_ 
2 

3 



1 

%/3 



(28) 



Of course, this other quartet of 77 and 77' decay constants 
can also be parameterized in terms of other two constants 
and two other mixing angles: 



f^^ = cos (/)NS ./l 



NS 



/,^ = -sin0s.fs , (29) 



f^,^ = sin 0NS /ns , ffj- = cos 4>s fs , (30) 
where /ns and /g arc given by the matrix elements 

{0\A^si^.)\ms{p))^^hsp''e-^P■- , (31) 



(0|Ag(a:)|r;s(p)) =7/sp"e-'^-" , (32) 

while {0\A^six)\rjs{p)) = = {0\A^{x)\ms{p)). 

In the NS-S basis, it is possible to recover a scheme 
with a single mixing angle (/) through the application of 
the Okubo-Zweig-Iizuka (OZI) rule [ii,!!!,!^^. For ex- 
ample, /ns/s sin((/)Ns — 4>s) differs from zero just by an 
OZI-suppressed term |28i. Neglecting this term thus im- 
plies (pNS = </>s- (Refs.^, 29,^] denote /ns, /s, 0ns, 0s 
by, respectively, fq, fs, 4>q, 4>s-) In general, neglecting the 
OZI-suppressed terms, i.e., application of the OZI rule, 
leads to the so-called FKS scheme [H, [H 113, which 
exploits a big practical difference between the (in prin- 
ciple equivalent) parameterizations ([23 |) -([24 |) and ([29]) - 
([5(7)) : while 9s and 6*0 differ a lot from each other and 
from the octet-singlet state mixing angle 9 « (6*8 -I- 6*0)72, 
the NS-S decay-constant mixing angles are very close to 
each other and both can be approximated by the state 
mixing angle: 0ns ~ 0s ~ 0- Therefore one can deal 
with only this one angle, 0, and express the physical ri-T]' 
decay constants as 



Jrj 
Jrj' 



Jrj' 



/ns cos 
/ns sin 



-/s sine 
/s cos ( 



1 

x/3 



2 
3 

1 

^/3 



(33) 

This relation is valid also in our approach, where 77 and 77' 
are the simple 77ns-?7s mixtures (|16|) . The FKS relations 
[2i,[2i[33 



/s - \l o/ns 



3^s 



= 0— arctan 



^/2/s 



/NS 



(34) 



A{^s(x) = -LA«^(x) + y^A°'^(x) 



1 



[u{x)-1^'-1^u{x) + J(x)7^75d(x)] , (26) 



/o - \/ o/ns + 



arctan 



V2/ 



NS 



(35) 
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equivalent to Eq. p3|) , were also shown [8| to hold in our 
DS approach. 

In our present DS approach, mesons arc pure qq 
BS solutions, without any gluonium admixtures, which 
are prominent possible sources of OZI violations. There- 
fore, our decay constants are calculated quantities, /ns = 
fuu = fdd = U and /s = fss, in agreement with the OZI 
rule. Our DS approach is thus naturally compatible with 
the FKS scheme, and we can use the j] and ?/ decay con- 
stants (j33[) with our calculated /ns = /ir and /g — fss in 
Shore's equations P^ - ipT]) . 



IV. RESULTS AND CONCLUSIONS 

All quantities appearing on the right-hand side of 
Eqs. ([TO)) - (PT|) . namely M.^, Mk, and Jk, are calcu- 
lated in our DS approach [I,[i3j[i3| (for the three dynam- 
ical models 0, [111, except the full QCD topological 
charge parameter A. Since it is at present unfortunately 
not yet known, we follow Shore and approximate it by 
the Yang-Mills topological susceptibihty xym- 

On the left-hand side of Eqs. ([T9l)-(I1T|), the model 
results for /ns = fir and /s = fss and Eq. ([33]) re- 
duce the unknown part of the four rj and rj' decay con- 
stants /", /",, /,^, and /^,, down to the mixing angle 
(j). The three Shore's equations (fTOl) - (PT|) can then be 
solved for </>, M„ and Af„/, providing us with the up- 



per three lines of Table |TT1 For each of the three dif- 
ferent dynamical models which we used in our previous 
DS studies 0, i, [13, [13, [13, HI, these results are dis- 
played for XYM = (175.7 MeV)4 as in Refs. [13, 



Inputs: 


from Ref. [8] 


from Ref. [10] 


from Ref. [13] 


1/4 


175.7 


191 


175.7 


191 


175.7 


191 


Mr, 


485.7 


499.8 


482.8 


496.7 


507.0 


526.2 




815.8 


931.4 


818.4 


934.9 


868.7 


983.2 


't' 


46.11° 


52.01° 


46.07° 


51.85° 


40.86° 


47.23° 


00 


1.84° 


7.74° 


1.39° 


7.17° 


-6.69° 


-0.33° 


Os 


-17.90° 


-12.00° 


-17.6° 


-11.85° 


-20.47° 


-14.11° 


fo 


108.8 


108.8 


107.9 


107.9 


101.8 


101.8 


fs 


122.6 


122.6 


121.1 


121.1 


110.7 


110.7 


rO 
J'l 


-3.5 


-14.7 


-2.6 


-13.5 


11.9 


0.6 


Jr,' 


108.8 


107.9 


107.9 


107.1 


101.1 


101.8 


rS 

Jv 


116.7 


119.9 


115.4 


118.5 


103.7 


107.4 




-37.7 


-25.5 


-37.6 


-24.9 


-38.7 


-27.0 



TABLE II: The results of the three DS models obtained 
through Shore's equations (|19p - (|2ip for the two values of xym 
approximating A: (175.7MeV)'' and (191MeV)*. Columns 1 
and 2: The results when the non-anomalous inputs for Eqs. 
((I9|-(l2l|, namely M.^,MK,/,r = /ns, fss = fs and /k, are 
taken from Ref. [3], which uses Jain-Munczek Ansatz interac- 
tion [^. Columns 3 and 4: The results for the non- anomalous 
inputs from Ref. [T^ using OPE-inspired interaction nonper- 
turbatively dressed by gluon condensates uM ■ Columns 5 and 
6: The results for the inputs from Ref. [l^] using the separa- 
ble Ansatz interaction [ij]- All masses and decay constants, 
as well as Xym' MeV, and angles are in degrees. 



for XYM = (191MeV)4 [p (adopted by Shore 
The lower part of the table, displaying various additional 
results, is then readily obtained through Eq. ([55]) and/or 
the useful relations which give /g, fo, ds^4' and 

6*0 - in terms of /ns = fw and fs = fss- Thus, unlike 
the mixing angles, /o and /§ do not result from solving 
of Eqs. (fT9|) -([21 ]) . but are the calculated predictions of a 
concrete dynamical DS model, independently of Shore's 
equations. 

For all three quite different (RGI [i,[ll| and non-RGI 
[3]) dynamical models which we used in our previous DS 
studies [3, [1, [13, [13, [13, [IHl , the situation with the results 
turns out to be rather similar. Similar results from vari- 
ous models mean that the usage of Shore's generalization 
in conjunction with the DS approach does not help one 
to discriminate between various dynamical models and 
so draw conclusions on the dynamics. This is not sur- 
prising, as it has been established [l|, [1, [3, [1, [3, [^ that 
while a successful reproduction of static properties and 
other low-energy meson phenomenology requires interac- 
tion modeling at low momenta, it is possible to achieve 
a satisfactory description of low-energy phenomenology 
for many forms of model interactions as long as their in- 
tegrated strength at low momenta {p^ < 1 GeV^) is suf- 
ficient to achieve a realistic DChSB. On the other hand, 
this (similarity of our results from the very different mod- 
els) has the advantage that our conclusions further below 
are not sensitive to the changes of the model dynamics. 

The most conspicuous feature of our results is that rj 
and rj' masses are both much too low when the weighted 
average xym = (175.7 ± 1.5 MeV)'' of Refs. [13, M, [11 
is used, in contrast to the results from the standard WV 
relation, displayed in Table [H If we single out just the 
highest of these values (191 MeV)^ [H), the masses im- 
prove somewhat. However, other results are spoiled - 
e.g., the mixing angle (j) becomes too high to enable agree- 
ment with the experimental results on rj, r/ — > 77 decays, 
which require (j) ^ 40° [l3l ■ 

When we turn to the lower parts of Tables |T] and [U 
where the results for the rj and rj' decay constants, and 
the corresponding two mixing angles 9q and 0g, are given, 
we notice a feature common to all our results, as well as 
Shore's (also given in Table |T|. The diagonal ones, /°, 
and /*, are all of the order of /^, being larger by some 
10% to 30%. The off-diagonal ones, /^, and /^, are, on 
the other hand, in general strongly suppressed. This is 
expected, as rj' is mostly singlet, and rj is mostly octet. 

To understand the dependence of the decay con- 
stants on the topological susceptibility xym (approximat- 
ing A), it is important to note that our /o, which in a 
bdl QCD bound-state calculation would be influenced by 
the gluon anomaly, presently is not, since it is calculated 
(same as /§) from the modeled meson qq substructure 
relying on RLA. In Tables I and II one therefore sees no 
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XYM-dependence of not only /§, but also of /o, since the 
difference between /o and /s is presently generated only 
by their different NS and S quark content. This feature 
is not only consistent with the FKS scheme, but is in 
fact a general characteristic of this scheme. Namely, due 
to the neglect of OZI- violating contributions [23|, in the 
SU(3) flavor symmetry limit one would have /ns = /s 
and /o = /s- (The DS approach will be able to obtain 
the gluon anomaly dependence of /o only when it man- 
ages to go beyond RLA, which is presently achieved only 
with schematic, very simplified (5-function-type interac- 
tions [13] ■) Thus, not only in the present DS calculation, 
but in fact in any application of the FKS scheme, the 
XYM-dependence of the four physical rj-rj' decay constants 
(l23ll - p4l) stems exclusively from the XYM-dependence of 
the mixing angles (f>, 9^ and Oq. Its origin, as explained in 
the previous section, is in the interplay of the anomaly 
with the flavor symmetry breaking. In fact, the FKS 
scheme is based on the assumption that the flavor sym- 
metry breaking is significantly more important than the 
OZI-violating contributions (arising beyond RLA in DS 
approach) . 

The feature that may be surprising is that Shore's 
results (which, to be sure, were obtained [l^, US] in quite 
a different way from ours) are more similar to our results 
obtained through the standard WV relation, than to our 
results obtained through Shore's Eqs. (fTOl) - (PT|l . 

To summarize: the present paper has explored a 
modification of the DS treatments of the V-v' complex 
employed in Kefs. 0,i,[l3,[ll- In Refs. 0,i,[ll,the 
value of the anomalous rjo mass shift was obtained by fit- 
ting, but Ref. [l^l improved the treatment by obtaining 
it from the lattice through the WV relation. Agener- 
alization of this relation was recently proposed [Ij, , 
and the purpose of the present paper is to tests the usage 
thereof in the bound-state, DS context, and compare the 
results with those from the standard WV relation. 

All in all, inspection and comparison of the results in 
Table HIl with the results (in Table |T| from the analogous 
calculations but using the standard WV relation to con- 
struct the complete 77-77' mass matrix, leads to the conclu- 
sion that the DS approach with the standard WV relation 
([TT)) is phenomenologically more successful, yielding the 
masses closer to the experimental ones. This may seem 



surprising, but one must be aware that we do not yet 
have at our disposal the full QCD topological charge pa- 
rameter A, and that we (along with Shore) had to use its 
lowest approximation, xym- This in general pre- 

cludes a consistently improved treatment in spite of 
the usage of Shore's relations. The problems with incon- 
sistencies in the 1/Nc counting may well cause spoiling of 
results, especially in an approach such as ours, where the 
?7 and 77' masses are not inputs, but predicted quantities. 
Our results thus add a new argument to the motivation 
for undertaking lattice calculations proposed by Shore 
[2^ 1 and aimed at proper finding the quantity A. Also, 
we should recall from Sections HI and HIl that the very usage 
of the RLA assumed that the anomaly is implemented on 
the level of the anomalous mass only, as a lowest order 
1/Nc correction 0, H, E [ll ■ Thus, with respect to 
the orders in l/Nc, using Shore's generalization in the 
present formulation of our DS approach may be less con- 
sistent than using the standard WV relation, which may 
well be the cause of its lesser phenomenological success. 

In spite of the lesser phenomenological success (than 
the standard WV relation) in the present context of 
bound-state DS calculations at zero temperature, the 
presently exposed usage of Shore's generalization will 
likely find its application at finite-temperature calcula- 
tions in the DS context. Namely, there it may help alle- 
viate the difficulties met due to the usage of the standard 
WV relation in the DS approach at T > 0, as discussed 
inRcf. [13. 
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